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Abstract. We present an XMM-Newton observation of the radio jet and diffuse halo of the nearby radio galaxy 
NGC 6251. The EPIC spectrum of the galaxy's halo is best-fitted by a thermal model with temperature kT ~ 1.6 
keV and sub-solar abundances. Interestingly, an additional hard X-ray component is required to fit the EPIC 
spectra of the halo above 3 keV, and is independently confirmed by an archival Chandra observation. However, 
its physical origin is not clear. Contribution from a population of undetected Low Mass X-ray Binaries seems 
unlikely. Instead, the hard X-ray component could be due to inverse Compton scattering of the CMB photons 
(IC/CMB) off relativistic electrons scattered throughout the halo of the galaxy, or non-thermal bremsstrahlung 
emission. The IC/CMB interpretation, together with limits on the diffuse radio emission, implies a very weak 
magnetic field, B « 1/iGauss, while a non-thermal bremsstrahlung origin implies the presence of a large number 
of very energetic electrons. We also detect X-ray emission from the outer (~ 3.5') jet, confirming previous ROSAT 
findings. Both the EPIC and ACIS spectra of the jet are best-fitted by a power law with photon index F ~ 1.2, 
fixed Galactic column density, and 1 keV flux F\ kcv = 2.1 njy. A thermal model is formally ruled out by the data. 
Assuming an origin of the X-rays from the jet via IC/CMB, as suggested by energetic arguments, and assuming 
equipartition implies a large Doppler factor (8 ~ 10). Alternatively, weaker beaming is possible for magnetic fields 
several orders of magnitude lower than the equipartition field. 
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1. Introduction 

NGC 6251 is a giant elliptical galaxy at redshift z=0.024. 
It is the host of a supermassive black hole with mass 
M BH ~ 4 - 8 x 1O 8 M (Ferrarese & Ford 1999), dy- 
namically estimated with HST. At radio wavelengths, the 
source is well studied. The radio source measures 1.2° 
across the sky with bright hotspots that lie toward the 
edges of its two bright radio lobes (Waggett et al. 1977). 
This structure is indicative of many powerful FanarofT & 
Riley (1974) type-II radio galaxies. However, its power at 
178 MHz (Waggett et al. 1977) is below the FanarofT & 
Riley division between FRI/II sources, so NGC 6251 ap- 
pears to be undcrluminous considering what is expected 
based on its morphology. A prominent radio jet appears 
to the N-W for 4.5', as well as a weak counterjet opposite 
of the nucleus (Perley et al. 1984). The radio jet appears 
narrow out to ~2', where the radio shows a bright knot 
and more diffuse emission (see Fig. 1). 
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Together with NGC 6252, which lies a few arcminutes 
to the North, NGC 6251 belongs to the outskirts of the 
cluster of galaxies Zw 1609.0+8212 (Young et al. 1979). 
As discussed in Birkinshaw & Worrall (1993), the cluster 
does not significantly affect the dynamical properties of 
NGC 6251. 

Previous X-ray imaging studies of NGC 6251 with 
ROSAT showed the presence of an unresolved nuclear 
source embedded in diffuse thermal emission associated 
with the galaxy's halo (Birkinshaw & Worrall 1993), 
with no contribution from the cluster gas. X-ray emission 
from the outer radio jet was also detected with ROSAT 
(Birkinshaw & Worrall 1993), with two compact knots at 
~ 4' and 6' from the core (Mack et al. 1997a). The latter 
authors favor a thermal origin for the X-rays on the ba- 
sis of energetic arguments. Optical emission from the jet 
region at ~ 20" was also claimed (Keel 1988), but never 
confirmed with HST. 

Comparison of the halo gas pressure to the internal 
jet pressure led Birkinshaw & Worrall (1993) to conclude 
that the jet cannot be confined by the atmosphere of the 
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galaxy. Furthermore, Mack et al. (1997a) proposed that 
the jet is magnetically confined. 

We observed NGC 6251 with XMM-Newton as part of 
a project aimed at studying the X-ray emission from the 
various components of the galaxy. With its better sen- 
sitivity and good resolution, the EPIC camera on-board 
XMM-Newton enables a detailed spectroscopic study of 
the nucleus and extended features in the 0.3-10 keV en- 
ergy range. The analysis of the nuclear properties is pre- 
sented in Gliozzi et al. (2003a; Paper I in the following). 
Here, we concentrate on the X-ray emission from the halo 
and the jet, using the EPIC data and archival Chandra 
and VLA data. In this paper, Ho = 75 km s _1 Mpc -1 and 
qo = 0.5 are adopted. With this choice, 1" corresponds to 
446 pc at the distance of NGC 6251 (97.6 Mpc). 

2. Observations and Data Reduction 

2.1. XMM-Newton 

XMM-Newton carries on-board three sets of CCDs in the 
focus of the telescope, the EPIC pn and the EPIC MOS1 
and MOS2 cameras. While the pn has a larger effective 
area and is thus optimized for higher signal-to-noise spec- 
troscopy, the MOS offers better imaging capabilities. Both 
the EPIC pn and MOS cameras operate in the energy 
range 0.3-10 keV. 

We observed NGC 6251 on 2001 March 26 for a dura- 
tion of ~ 41 ks with the EPIC pn, and for ~ 49 ks with 
EPIC MOS1 and MOS2. The EPIC pn camera was oper- 
ated with the thin filter, and the MOS cameras with the 
medium filter (to avoid pileup of the bright nucleus; Paper 
I) . Details of the reduction procedure were given in Paper 
I. In summary, the data were reprocessed and screened 
with the XMM-Newton Science Analysis Software (SAS 
v.5.3.3) to remove known hot pixels and other data flagged 
as bad; only data with FLAG=0 were used. Investigation of 
the full-field light curves revealed the presence of a pe- 
riod of background flaring at the end of the observation. 
These events were screened, reducing the effective total 
exposures to ~ 36 ks for the EPIC pn and ~ 43 ks for 
the MOS cameras. Background data were extracted from 
source-free circular regions on the same chips containing 
the source. There were no signs of pile-up in the pn or 
MOS cameras according to the SAS task epatplot. 

The X-ray spectrum of the jet was extracted from an 
elliptical region with semi-axes of 40"and 12", located at 
~ 230" from the nucleus. The total count rate of the jet in 
the above extraction region in the energy range 0.5-9 keV 
is (9.3 ± 1.4) x 10~ 3 c/s with EPIC pn, (2.8 ± 0.5) x 10~ 3 
c/s with MOS1, and (1.5 ± 0.5) x 10~ 3 c/s with MOS2. 

Two different regions were used to investigate the spec- 
tral properties of the halo: 1) an annular region with radii 
35-215"; 2) a series of annular regions, with an incre- 
ment of the radius of 30", with inner radii ranging from 
35" to 215". In both cases, the radio jet and the point 
sources visible in the EPIC MOS image were removed. 
Conservatively, a region of the same size as the jet but at 
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Fig. 1. (a, Top) EPIC pn image of NGC 6251 in 0.3- 
10 keV and (b, Bottom) Chandra ACIS-I image in 0.3-8 
keV. The resolution of the EPIC image is 10", while the 
resolution of the ACIS image is 2". In the two panels, the 
1.4 GHz VLA data are overlaid. X-ray emission from the 
nucleus, diffuse halo, and outer jet is apparent. 



the opposite azimuth (where a possible counterjet could 
be present at X-rays) was also excised. The count rate of 
the halo in an extraction annulus with radii 35-215" in 
the energy range 0.5-2 keV is 0.37 ± 0.01 c/s with EPIC 
pn, 0.102 ± 0.003 c/s with MOS1, and 0.110 ± 0.003 c/s 
with MOS2. 

Spectral analysis was performed using XSPEC v.ll. The 
EPIC pn and MOS spectra were fitted jointly, leaving the 
relative normalizations free to vary to account for inter- 
calibration effects. The data were rebinned with a mini- 
mum of 20 counts per bin, in order to apply the \ 2 statis- 
tics, and fitted using the latest calibration files provided by 
the XMM-Newton Guest Observer Facility. Improvement 
of the fit when free parameters were added to the fit was 
evaluated using the F-test, assuming as a threshold for 
improvement Pf=95%. 
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2.2. Chandra 

A 25.4 ks archival Chandra image of NGC 6251 was inves- 
tigated to derive better spatial constraints on the jet and 
halo. Chandra observed NGC 6251 on Sept 11, 2000 for 
25.4 ks with ACIS-I (PI Kerp). In order to study the X- 
ray emission from the jet, the source was positioned with 
the nucleus on the gap between chips II and 12. While 
this strategy prevents a detailed analysis of the nuclear 
properties, the Chandra data in principle can be useful to 
investigate the spatial properties of the X-ray halo and 
jet. 

We reduced the ACIS data using standard procedures 
with the CIA02.3 software. We restricted our analysis to 
the 0.3-8.0 keV energy range, where the instrument is 
better calibrated and the background is negligible. We in- 
spected the light curves of the background to search for 
possible background fluctuations. No flares have been de- 
tected. Before extracting the radial profiles of the sources, 
we removed field sources from the images (see Appendix) . 
To find the field sources we used the tool wavdetect with 
the default value (10~ 9 ) for the threshold for identifying 
a pixel as belonging to a source. 

Using the tools dmextract and dmtcalc, we extracted 
the radial profiles on a scale of 4'. The jet and halo spec- 
tra were extracted in the same regions as for EPIC (see 
above), removing the point sources and, in the case of 
the halo, the jet and possible counterjet regions. In order 
to evaluate the instrumental response to a point source, 
we used the method explained in Donato et al. (2003). 
Although a script that automatically performs the same 
task (ChaRT) has been released by the Chandra X-ray 
Center, this script failed in the creation of the instrumen- 
tal PSF, since the source is positioned in the gap between 
chips (photons are created but not detected). 

2.3. Radio 

In order to explore the spatial coincidence of X-ray fea- 
tures with the known radio jet, we obtained data from 
the archive of the NRAO 1 Very Large Array (Thompson 
et al. 1980). We chose a L-band observation from August 
15, 1995 which used the full VLA to collect over 3 hours of 
data "on-source" in its highest resolution A-configuration 
in spcctral-linc mode. This data has been published by its 
original investigators (Werner et al. 2001; Werner 2002) 
and was processed similarly. In brief, we used AIPS (Bridle 
& Greisen 1994) for the initial calibration utilizing scans 
of 3C 286 to set the flux densities to the VLA 1999.2 
scale. The nearby source 1800+784 was used as the phase 
calibrator for NGC 6251, and bandpass calibration per- 
formed using 3C 286. The inner 22/31 channels were av- 
eraged to produce a single —9 MHz bandwidth contin- 
uum dataset centered at about 1.4 GHz which was then 
self-calibrated in the Caltech DIFMAP package (Shepherd 

1 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 



et al. 1994). The final images were corrected for primary 
beam response using the AIPS task PBCOR. This correc- 
tion amounts to up to an —5% correction near the ~4' 
distant X-ray feature in the jet. 

3. Results 

The EPIC pn image of NGC 6251 in 0.3-10 keV is shown 
in Fig. la, with the VLA contours at 1.4 GHz superposed. 
The image was binned by a factor 32 and smoothed with 
a Gaussian function with width a=6". The final resolu- 
tion is -10". The smoothed ACIS-I image in 0.3-8 keV 
is shown in Fig. lb. This image was binned by a factor 
4 and smoothed with the CSM00TH with final resolution 
2". The radio images overlaid in the figures were restored 
with beamsizes matching the resolutions of the respective 
X-ray images. 

As apparent in Fig. la-b, various components con- 
tribute to the total X-ray emission of the galaxy: the nu- 
cleus, a diffuse isotropic halo associated with the elliptical 
galaxy (e.g., Birkinshaw & Worrall 1993), and jet emission 
at — 4' from the core. Also apparent in the Chandra im- 
age are several point sources scattered as far as 4-5' from 
the nucleus (see Appendix). Using the Chandra data, we 
determined that their contribution to the total X-ray emis- 
sion in the EPIC extraction radius of the halo is ~ 5%. 

3.1. The X-ray halo 

To investigate the spatial properties of the halo, we ex- 
tracted brightness radial profiles from EPIC MOS2 on- 
board XMM-Newton and ACIS-I on-board Chandra. The 
choice of the MOS2 detector is motivated by geometrical 
considerations, namely, the fact that the source is centered 
on its CCD1, and its higher sensitivity than MOS1, offer- 
ing the best chance to study the extended emission. The 
EPIC radial profile in the total energy band 0.3-10 keV 
is shown in Fig. 2a, together with the instrumental PSF 
and the background, while the ACIS-I profile is shown in 
Fig. 2b. The MOS2 profile is fully consistent with the in- 
strumental PSF (the MOS2 PSF was obtained using the 
analytical formula from Ghizzardi 2001) plus background, 
whereas the ACIS profile suggests the presence of an ex- 
cess between 20" and 100". However, the limited signal-to- 
noise ratio of the ACIS data hampers a more quantitative 
spatial analysis. We fitted the Chandra profile with a (3 
model; the fitted parameters are unconstrained. 

Fig. 2c shows the energy-dependent EPIC profiles at 
soft (0.3-2 keV) and hard (2-10 keV) X-rays. The soft 
and hard X-ray profiles are similar, suggesting that the 
the hard X-ray emission originates within the galaxy. 

EPIC pn and MOS spectra for the halo were ex- 
tracted as described above. We first fitted the total in- 
tegrated spectrum with a thermal plasma model (apec in 
XSPEC), leaving as free parameters the temperature, fcT, 
the abundance, Z (in solar units), and the normalization. 
All components are absorbed by a column density fixed 
to the Galactic value, N§ al = 5.5 x 10 20 cm" 2 (Dickey & 
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Lockman 1990). The fit gives \r = 1-2/2913. Inspection 
of the residuals shows the presence of excess flux above 
the model at energies ~ 3 keV (Fig. 3a). Adding a power 
law to the model gives x? =1-10/29 11, improved at Pp > 
99% with respect to the fit with the thermal alone. When 
the hard tail is modeled with a thermal brcmsstrahlung 
component, the fit is worse (x 2 = 1.14/2911) than with a 
power law. The best-fit temperature is 199 keV, effectively 
mimicking a flat power law. 

The fitted parameters of the thermal plus power law 
model are: kT = 1.57 ± 0.09 keV, Z = 0.10 ± 0.03Z Q , 
and photon index T = — 0.08to^ 8 - The confidence con- 
tours for the temperature and abundance of the thermal 
component are shown in Fig. 3b. The observed 0.5-2 keV 
flux of the thermal is Fo.5-2 koV ~ 3.1 - 4.7 x 10~ 12 erg 
cm~ 2 s , and the intrinsic (absorption-corrected) 0.5-2 
keV luminosity is Lo.5-2 keV ~ 4 - 6 x 10 42 erg s _1 . The 
power law accounts for ~ 20% of the total X-ray emission 
in 0.5-2 keV. Its observed flux and intrinsic luminosity in 
2-10 keV are F 2 _i kc v ~ 1-7 - 2.6 x lO" 11 erg cm~ 2 s"\ 
and L 2 -io kcv ~ 1-8 — 2.7 x 10 42 erg s _1 , respectively. 

We investigated whether the hard power law is related 
to the contaminating wings of the PSF (the nucleus is 
a bright X-ray emitter; Paper I). We extracted the halo 
spectra in annuli of increasing radii, as described in Sect. 
2.1. If these spectra are fitted with a thermal model only, 
the fit is acceptable (x 2 = 0.9 — 1.2) but the fitted tem- 
perature is very high, kT = 4 — 5 keV, much larger than 
typically found for the halo of radio galaxies, kT ~ 1 
keV (e.g., Worrall, et al. 2001). Inspection of the residuals 
shows that a power law is present even in the spectrum ex- 
tracted in the outermost annulus between 155" and 215". 
If the power law is added, the fits improve with high sig- 
nificance, Pp ~ 99%; its photon index is not well defined 
and ranges between -0.7 and 1.2. Within the uncertain- 
ties, there is no variation of either the temperature or the 
abundance of the thermal component with distance from 
core. The temperature varies in the range 1.3-1.9 keV and 
the abundance in 0.05-0.2 (90% uncertainties). Moreover, 
the nucleus spectrum has a substantially steeper slope, 
T ~ 1.9 (Paper I), than the hard component. It is thus 
unlikely that the power law represents the contribution of 
the PSF wings. 

A non-thermal component is also required to fit the 
Chandra spectrum of the halo. The latter was extracted 
in a similar region as for EPIC, an annulus of radii 35" and 
215", after removing the jet region and the point sources. 
Spectra were extracted for each CCD and fitted jointly. 
A fit with a thermal only is good (x r 2 = 0.96/209), but 
yields an unacceptably high temperature kT ~ 10 keV for 
Z = 0.08. Adding a power law causes the x 2 to decrease 
by A% 2 =35, significant at > 99.9% confidence. The fitted 
temperature then converges towards the EPIC value with 
kT = 1.25±g;|? keV. 

To investigate the spatial distribution of the hard X- 
ray component, we extracted soft (0.5-2 keV) and hard (2- 
8 keV) ACIS images. The hard X-ray image is dominated 



by the nucleus point source and there are not enough 
counts in the diffuse emission to study its spatial prop- 
erties. 

In summary, the EPIC spectrum of the diffuse halo in 
NGC 6251 is best-fitted by a thermal model with temper- 
ature ~ 1.6 keV and sub-solar abundances. There are no 
gradients of cither parameter with distance from the core. 
The data also require a hard X-ray component above 3 
keV, which is best-fitted by a power law with a flat slope, 
albeit poorly determined. The hard component, which ac- 
counts for 20% of the total 0.5-10 keV luminosity of the 
halo, is also detected in the Chandra spectrum. 



3.2. X-ray emission from the jet 

Fig. la-b shows that weak X-ray emission is detected at 
~ 4' from the core in both the EPIC and ACIS images, 
coincident with the bright radio jet. To assess the signifi- 
cance of the jet detection we compared the average counts 
found in the elliptical jet region with the counts found 
in six identical elliptical regions randomly located around 
the jet. The significance of the jet detection is of the or- 
der of 10cr. We thus confirm the previous claim based on 
ROSAT of X-ray emission from the outer jet (Mack et al. 
1997a). The X-ray emission of the jet is diffuse in both 
the EPIC and ACIS images, with no compact knots at 
the resolution of ACIS-I. 

We extracted X-ray spectra of the jet with EPIC. The 
pn and MOS spectra were fitted with a single power law 
model with Galactic absorption, 5.5 x 10 20 cm~ 2 . This fit is 
acceptable, x 2 =0.98 for 102 de grees of freedom, and gives 
a best- fit photon index T = 1.15+oJg. An equally accept- 
able fit was obtained with a thermal model; however, the 
fitted temperature is kT = 63 keV and the abundances 
near zero, effectively mimicking a power law. A fit with a 
thermal plus a power-law model does not lead to signifi- 
cant improvement. Using the best-fit power law model, the 
observed 2-10 keV flux is F 2 - W kc v = (1-6 - 3.1) x 10~ 14 
erg cm~ 2 s _1 , and the 1 keV flux is F\ koV = 2. 31^6 nJy. 

A similar result is obtained with Chandra. The ACIS-I 
spectrum of the jet was extracted in a 38" x 22" elliptical 
region. In the energy range 0.3-8 keV, about 109 counts 
were derived, allowing a crude spectral analysis. Again, 
the model that best-fits the ACIS data is a power law, with 
r = 1.06ig;38 and 1 kcV flux F i kcv = 3.ll?;o nJy. These 
parameters are consistent with XMM-Newton within the 
uncertainties. As for the EPIC spectra, a thermal model 
is again formally acceptable but the fit converges to a 
temperature kT ~ 68 keV and near-zero abundances. 

In summary, X-ray emission from the jet at 4' from the 
core is detected with both XMM-Newton and Chandra. 
The X-ray spectrum of the jet is better fitted by a flat 
(r ~ 1.2) power law. 
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4. Discussion 

4.1. The X-ray halo 

NGC 6251 is embedded in a diffuse X-ray halo extend- 
ing ~ 100", or ~ 45 kpc from the core. The fitted 
temperature of the halo, kT ~ 1.6 kcV, is significantly 
larger than derived from ROSAT PSPC data, - 0.5 keV 
(Birkinshaw & Worrall 1993). This is not surprising, as 
NGC 6251 is thought to host a cooling flow in the nu- 
clear region (Birkinshaw & Worrall 1993), which was ex- 
cluded from our extraction region. Indeed, spectral fits 
to the EPIC spectrum of the nucleus requires a thermal 
component with a temperature kT ~ 0.55 (Paper I), in 
agreement with the ROSAT data. 

A new result of our analysis is the detection of a hard 
X-ray component in the spectrum of the halo above 3 keV. 
The hard tail is best described by a power law, contribut- 
ing ~ 20% of the total 0.5-10 keV X-ray flux. A thermal 
model for the X-ray tail is inconsistent with the data. The 
hard non-thermal component is independently confirmed 
by an archival Chandra observation. 

Inspection of the Chandra image (Fig. lb) shows that 
the hard component in the EPIC spectrum cannot be due 
to the collective contribution of the detected off-axis X-ray 
point sources. The hard X-ray component is still present 
in the ACIS data when the point sources are excised. The 
hard X-ray profile follows closely the soft profile (Fig. 2c), 
indicating that the hard X-ray emission is related to the 
galaxy halo. 

We thus investigated the possibility that the hard X- 
ray component could be due to the emission of a popu- 
lation of undetected Low-Mass X-ray Binaries (LMXBs). 
Based on recent Chandra observations of nearby ellipticals 
(e.g., NGC 1316, Kim & Fabbiano 2003; NGC 5128, Kraft 
et al. 2003; NGC 1399, Angelini et al. 2001), this contribu- 
tion could be substantial. First, we examined the nature 
of the detected off-axis point sources in Table 1 . Using the 
count rates reported in the Table, and assuming a power 
law model with T = 1.5 — 2.0 and Galactic absorption, we 
estimate that the threshold X-ray flux in 0.3-8 keV for 
detecting point sources is 3 x 10~ 15 erg cm -2 s _1 . This 
corresponds to a luminosity L .3-8 kcV ~ 3 x 10 39 erg s _1 , 
at least one order of magnitude larger than for LMXBs. 

Based on the luminosity functions of field sources from 
X-ray deep surveys (Kim et al. 2003), in a 4' radius cen- 
tered on NGC 6251 we expect 14 ± 4 serendipitous sources 
at soft X-rays, consistent with our measured rate of 10 
sources (Table 1). It is thus likely that the off-axis point 
sources in Table 1 are field galaxies, although we cannot 
exclude the possibility that some are Ultra-Luminous X- 
ray sources, which indeed were found in other ellipticals 
(e.g., Jeltema et al. 2003). 

If a population of LMXBs is indeed present in 
NGC 6251, it is undetected in our XMM-Newton and 
Chandra images. However, its presence could affect the 
integrated spectrum. To check for this, we added a new 
component to the best-fit model in Sect. 3.1 mimicking 



the integrated contribution of the LMXBs. Following Kim 
& Fabbiano (2003), we modeled this component with ei- 
ther a steep (r ~ 1.8) power law, or a kT ~ 5 keV thermal 
model. The addition of either the power law or thermal 
does not improve the fit to the halo spectrum. Moreover, 
the spectral parameters of the best-fit model components 
(Sect. 3.1) are basically unaffected. 

An independent indicator of the relative contribution 
of LXMBs to the X-ray emission of the galaxy is given 
by the ratio of the X-ray-to-optical luminosities, Lx/Lb, 
where Lb is in solar units. Kim et al. (1992) showed that 
the LMXB contribution decreases with increasing Lx/Lb- 
For NGC 6251, we derive log(L x /L B ) =±32.3, which ac- 
cording to Kim et al. (1992) indicates that the LMXB con- 
tribution to the hard X-ray emission is negligible. This is 
confirmed by the later analysis of Matsumoto et al. (1997), 
who analyzed ASCA observations of 12 early-type galax- 
ies. The X-ray halo luminosity of NGC 6251 is two orders 
of magnitude larger than predicted by the Lx — Lb rela- 
tionship in Fig. 6 of Matsumoto et al. (1997). We assumed 
the halo luminosity given in Sect. 3.1, and an optical mag- 
nitude B = 13.64 mag (de Vaucouleurs et al. 1991). Based 
on this evidence, we conclude that it is unlikely that the 
hard X-ray tail in the EPIC spectrum of the NGC 6251 
halo is due mostly to unresolved LMXBs. 

The low abundances we detect for the halo are puz- 
zling. While in some sources (e.g., NGC 1316; Kim & 
Fabbiano 2003) they seem to be related to LMXB contam- 
ination of the integrated diffuse emission, this is not the 
case for NGC 6251. Deeper X-ray observations designed 
to probe the structure of the diffuse thermal emission are 
needed. 

We explored alternative origins for the hard X-ray 
component. There are theoretical reasons to believe that 
non-thermal halos should be present in clusters/groups 
containing radio galaxies, as the latter are important 
sources of energetic particles (e.g., Ensslin et al. 1997). 
The energy stored in intcrgalactic magnetic fields can be 
large and will last long after the activity in the nucleus 
is turned off (Kronberg et al. 2001), heating the ICM 
through dissipation processes. Interestingly, in NGC 6251 
the ~ 35" halo temperature, kT ~ 1.6 keV, is hotter 
than in other "normal" giant elliptical galaxies (e.g., Xu 
et al. 2002) and to other less powerful radio galaxies (e.g., 
NGC 4261; Gliozzi et al. 2003b), where kT ~ 0.5 - 0.7 
keV. This lends credit to the idea that the AGN can in- 
ject particles and energy heating the ISM. 

The hard component detected in NGC 6251 is sim- 
ilar to the one previously observed with ASCA in the 
nearby group of galaxies HCG 62 (Fukazawa et al. 2001). 
The latter authors find that either a power law with pho- 
ton index T = 0.8 - 2.7, or a very hot (kT > 6 keV) 
bremsstrahlung best-fits the data at energies > 4 kcV, 
in both cases accounting for ~ 20% of the total 0.5-10 
keV flux of the halo. Fukazawa et al. (2001) conclude that 
the most likely origin of the hard X-ray component is ei- 
ther inverse Compton scattering of the Cosmic Microwave 
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Background photons (IC/CMB) on a population of low- 
energy (7 ~ 10 3 ) electrons, if the magnetic field of the 
ICM is less than 1/iGauss, or non-thermal bremsstrahlung 
from a population of energetically important relativistic 
electrons. More recently, non-thermal diffuse emission was 
detected in a deep Chandra image of the distant (2=3.8) 
radio galaxy 4C 41.17 (Scharf et al. 2003), and interpreted 
as IC on the local CMB and FIR radiation fields. Following 
Fukazawa et al. (2001), a possible origin for the X-ray tail 
in NGC 6251 is IC/CMB off a population of old electrons. 
This possibility is particularly attractive considering that 
NGC 6251 has an active radio core and a long jet which 
disrupts at ~ 4' from the core, where diffuse X-ray emis- 
sion is detected. If the jet magnetic field is disrupted and 
becomes unable to confine the electrons, the latter could 
have escaped and diffused in the ICM of the galaxy. The 
relativistic electrons would then lose energy via IC scat- 
tering of the CMB photons, whose energy density at the 
redshift of NGC 6251 is non-negligible (see below). 

However, the very flat slope of the hard X-ray compo- 
nent we measure from both XMM-Newton and Chandra 
spectra, Y ~ (albeit with large uncertainties), posits a 
problem for the IC/CMB interpretation. In fact, this im- 
plies that the minimum Lorentz factor of the scattering 
electrons is large, 7 m i n ~ 3000. We would thus be observ- 
ing the low-energy Compton tail emission (characterized 
by a spectral slope a = — 1) from a very energetic electron 
population. 

Kronberg et al. estimate for NGC 6251 a minimum 
total energy content of the X-ray halo of ~ 4.2 x 10 59 ergs 
and a minimum halo magnetic field of ~ 4/iGauss. An 
immediate consequence of such large magnetic field and 
energetic electrons is that one should observe large radio 
fluxes from the halo. 

A deep 326 MHz Westerbork Synthesis Radio 
Telescope (WSRT) map (Mack et al. 1997b) shows no ob- 
vious sign of such a radio halo surrounding this galaxy. 
The only diffuse emission in the image is detected in the 
actual radio lobes near the edges of the source. We esti- 
mate an integrated flux limit of 275 mJy (3 sigma) for 
the halo out to 215" (the extent of the XMM-Newton ex- 
traction aperture) from measuring the RMS in this re- 
gion of the digital FITS version 326 MHz map of the 
Mack et al. (1997b) which was available to download from 
DRAGN web-site 2 . Fig. 4 shows the plot of the predicted 
synchrotron flux at 326 MHz versus the halo magnetic 
field for two different values of the radio spectral index. 
In this calculation we assume that the observed hard X- 
ray power-law is due to the low-energy tail of electrons 
at 7 m in ~ 3000. Above an energy E ~ 10 keV the (un- 
observed) emission would be a power-law with slope a, 
produced by electrons with 7 > 7 m ; n . In order to repro- 
duce the observed radio flux (horizontal dotted line), a 

magnetic field -Bhaio ~ 10~ 3 /iGauss is required. This is a 
factor 10 or more smaller than generally measured for ra- 

2 J. P. Leahy, A. H. B ridle, & R. G. Strom, 
:/ /www .jb.man.ac.uk/atlas/ 



dio galaxies (Feretti et al. 1995). However, measurements 
of rich clusters with BeppoSAX and RXTE suggest that in- 
tergalactic magnetic fields can be very weak (e.g., Valinia 
et al. 1999). Note that, however, for magnetic fields be- 
low B ~ 1/xG, if electrons have a high energy cut-off at 
7max < 10 4 the maximum frequency of the synchrotron 
emission will be located below ~ SOOBi^cTmax 4 MHz. 
Therefore if the high-energy non-thermal electron popula- 
tion has a narrow electron distribution, 3000 < 7 < 10 4 , 
the synchrotron emission would be unobservable at ~ 300 
MHz even with moderately low magnetic fields, B ~ 1/xG. 
A similar solution (IC/CMB scattering off low-energy elec- 
trons, emitting synchrotron radiation well below the ob- 
served frequencies) was independently proposed for the 
high-redshift (z = 4.3) quasar GB1508+5714, where an X- 
ray halo with no radio counterpart was recently detected 
(Yuan et al. 2003). 

Alternatively, as discussed by Fukazawa et al. (2001) 
for HCG 62, the X-ray tail could be due to non-thermal 
bremsstrahlung involving subrelativistic but suprather- 
mal electrons (e.g., Kempner & Sarazin 2000). Assuming 
the ratio of thermal to non-thermal luminosity of 0.2, 
Fukazawa et al. derive that the energy density of non- 
thermal electrons would be ~ 0.6 times that of ther- 
mal electrons. Similar conclusions apply in the case of 
NGC 6251, where the non-thermal component accounts 
for a similar fraction of the diffuse X-ray flux. Kempner 
& Sarazin (2000) discuss a model where the non-thermal 
X-ray tail in the cluster Abell 2199 is produced by a popu- 
lation of electrons distributed in energy as a power law of 
index /i. For a steep distribution (fi ~ 3.5), the X-ray tail 
can be approximated by a power law with Y ~ 1 + 0.5/x. In 
our case, we would derive /1 ~ — 2, implying the presence 
of a large number of very energetic electrons. As these elec- 
trons would lose energy via synchrotron, again one would 
expect to observe diffuse radio emission. 

In summary, we find evidence for the presence of a 
non-thermal high-energy tail in the halo of NGC 6251 
with both XMM-Newton and Chandra. The physical ori- 
gin of this component, however, is not clear. A possible 
explanation - IC/CMB off electrons in the halo - implies 
very weak magnetic fields. Non-thermal bremmstrahlung 
requires the presence of a large number of very energetic 
electrons. 

4.2. The X-ray jet 

The EPIC and ACIS data confirm the previous findings 
from ROSAT that the outer jet emits X-rays. The X- 
ray spectrum of the jet is best described by a power law 
with r ~ 1.2, while a thermal model is formally ruled 
out by the data. Fig. 5 shows the radio-to-X-ray spec- 
tral energy distribution of the jet. The radio flux at 1.4 
GHz was measured (see above) using the archival VLA im- 
age and an extraction region similar to the EPIC extrac- 
tion region. The radio spectrum between 1.4 and 5 GHz 
throughout the jet was determined to be ct r = 0.64 ± 0.05 
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by Perley et al. (1984). We adopted this spectrum and 
used our 1.4 GHz measurement to calculate representa- 
tive flux values at 1.66 and 4.885 GHz for the region of 
interest. The open bowtie represents the best-fit EPIC 
spectrum. Extrapolating the radio spectrum to the X-ray 
region (dotted line in Fig. 5) overestimates the X-ray flux 
by more than 2 orders of magnitude. Moreover, the X-ray 
spectrum is much flatter than the lower-energy extrapola- 
tion. This suggests that the X-ray emission is not the high- 
energy tail of the synchrotron responsible for the longer 
wavelengths, and that a different spectral component is 
needed to account for the X-rays. 

Based on the ROSAT PSPC detection, Mack et al. 
(1997a) argue in favor of a thermal origin for the X-ray 
emission. However, a thermal origin for the X-ray emission 
is ruled out independently by the EPIC and ACIS spectra 
(see Sect. 3.2). We also note that there is no evidence 
for internal depolarization from the radio (Perley et al. 
1984). This makes a thermal origin for the X-rays unlikely. 
However, vis-a-vis the limited signal-to-noise ratio of the 
EPIC /ACIS spectra, the possibility that at least a fraction 
of the X-rays are due to thermal emission from ambient 
gas remains open. 

A likely possibility is inverse Compton (IC) scattering. 
Seed photons for IC can be provided by 1) the synchrotron 
photons themselves (SSC process); 2) the beamed emis- 
sion from the inner jet; 3) the isotropic emission of the 
AGN; and 4) the Cosmic Microwave Background radia- 
tion (CMB). At a projected distance of 94 kpc from the 
nucleus, the Narrow Line Region and dusty torus are not 
a significant source of photons. 

We compared the radiation densities (in the blob co- 
moving frame) at the blob location for processes 1-4 in 
order to evaluate the dominant source of photons. The 
SSC luminosity was obtained integrating the radio spec- 
trum in Fig. 5 up to 10 12 Hz, L s ~ 1.9 x 10 41 erg s" 1 . 
Assuming a size for the X-ray emission region of 10" 
(the EPIC resolution), the synchrotron radiation density 
is u s ~ 2.8 x 10~ 15 (5~ 4 erg cm~ 3 . Here 5 is the jet Dopplcr 
factor, defined as S — [Fl(1 — /3cos^)] _1 , with Tl the jet 
Lorentz factor, fi the plasma speed, and 9 the jet inclina- 
tion to the line of sight. 

The contribution from the inner jet (e.g., Celotti et al. 
2001) can be evaluated as u in ~ 3 x 10~ 15 L in , 46 , Tr^r^ 2 
erg cm -3 , where L ln is the typical observed luminosity of 
blazars (e.g., Ghisellini et al. 1998) and Ti n is the Lorentz 
factor of the inner jet. 

For the AGN luminosity, we used the bolometric lu- 
minosity from Paper I, Lagn ~4x 10 43 erg s _1 . At the 
location of the X-ray jet, the radiation density due to the 
isotropic radiation of the AGN is magn ~ f -3 x 1Q~ 15 T^ 2 
erg cm~ 3 . 

The CMB radiation density scales like (1 + z) 4 and r£. 
We find mcmb ~ 4.4 x lO -1 ^ erg cm~ 3 . This is at least 2 
orders of magnitude larger than u s , u- mi and «agn- Thus, 
the CMB is the dominant source of seed photons for the 
IC process. Note that assuming significant beaming in- 



creases the importance of wqmb and decreases that of the 
other processes. To examine this in more detail, we plot in 
Fig. 6 the relationship between the Doppler factor, S, and 
the magnetic field, B, for the IC/CMB process (dashed 
line), the SSC process (continuous line), and the equipar- 
tition condition (dotted line). Assuming equipartition, the 
IC/CMB process accounts for the X-ray emission of the 
outer jet if 5 ~ 10 and B ~ B cq ~ 2 /xGauss. Unless the 
magnetic field is extremely small (<< 0.1 /xGauss), SSC 
requires de-beaming. 

Assuming the equipartition magnetic field, the ra- 
diative time of the synchrotron electrons at 1.4 GHz is 
t r ~ 10 s yr. This is much larger than the light crossing 
time of the radio/X-ray emission region, t croS s ~ 2.4 x 10 4 
yr, accounting for the diffuse morphology at radio and 
X-rays. 

The assumption of IC/CMB in equipartition requires 
a strongly beamed (5 ~ 10) jet at 3.5' (so at least 93.7 
kpc) from the core. For a jet inclination ~ 30°, as 
expected from unifying schemes (e.g., Urry & Padovani 
1995), this implies large Lorentz factors for the emitting 
plasma. Alternatively, weaker beaming (S << 10) is al- 
lowed for a magnetic field at least one order of magnitude 
away from equipartition (Fig. 6). 

We thus explored other possibilities accounting for 
the concave radio-to- X-rays spectrum in jets, such as the 
model proposed by Dermer & Atoyan (2002). The latter 
authors interpret the continuum from radio to X-rays as 
synchrotron emission from a single electron distribution, 
whose shape, due to the different IC cooling rate suffered 
by electrons with different energies, produces the optical 
"valley". The model assumes that the electron cooling is 
dominated by IC scattering of the CMB photons. Low- 
energy electrons scatter the CMB in the Thomson limit, 
while the scattering by high-energy electrons, with Lorentz 
factor above 7kn ~ toc2 /^cmB' w * u be suppressed by 
the Klein-Nishina decline of the cross section. Due to the 
low cooling rate of the high-energy electrons, the result- 
ing electron energy distribution will be characterized by a 
flattening, at large Lorentz factors, which in turn produces 
a "bump" in the vF(y) synchrotron spectrum. 

We have calculated the expected spectrum implement- 
ing the simple analytical formulas provided by Dermer 
& Atoyan (2002). We have assumed that the emission is 
very weakly beamed (6 = 1). The result is presented in 
Fig. 7. Note that, due to the assumed dominance of IC 
losses, the magnetic field is constrained to be less than 
B ~ (8ttUcmb) ' 5 ~ 10~ 6 G (assuming no beaming). In 
our case this implies that the emitting plasma is strongly 
far from equipartition, since the number of electrons re- 
quired to produce the observed flux implies U e ~ 10 -9 erg 
cm~ 3 (weakly dependent on the minimum Lorentz factor 

7min)- 

To conclude, we cannot identify a simple viable way to 
account for the observed X-ray emission from the jet, since 
all the possibilities explored above present important dif- 
ficulties. However, it seems possible to conclude that the 
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jet is not confined by the external gas. In fact, assuming Appendix: X-ray point sources in the Chandra field 
a halo size of 200" or 2.86 x 10 23 cm, the density of the 
gas is n e ~ 2 x 10~ 3 cm~ 3 , implying a gas pressure of 



-Pext ~ 5 x 10~ 12 erg cm~ 3 . If the plasma in the jet is not 
too far from equipartition with the magnetic field (see dot- 
ted line in Fig. 6), this value is smaller than or very close 
to the pressure inside the jet. On the other hand, even if 
the magnetic field is well below equipartition, the pres- 
sure of the relativistic electrons necessary to account for 
the observed X-ray emission, both in the IC/CMB model 
and in the Dermer & Atoyan model, is much larger then 
the pressure in the external gas. We conclude the X-ray 
jet is overpressured with respect to the external medium, 
and cannot be confined by the halo, in agreement with 
previous findings based on ROSAT (Mack et al. 1997a, 
Birkinshaw & Worrall 1993). 



5. Summary 

We presented a 40 ks XMM-Newton EPIC observation 
of the giant elliptical NGC 6251, host of an FRI/II radio 
galaxy. An archival 25 ks Chandra exposure was also used. 
Our results can be summarized as follows: 

— X-ray emission from the region of the 4' jet is detected, 
as well as diffuse emission from the galaxy's halo. 

— The X-ray spectrum of the halo is described by a ther- 
mal model with temperature kT ~ 1.6 keV and sub- 
solar abundances, Z ~ 0AZ & . 

— A hard X-ray component is required to fit the EPIC 
spectrum of the halo above 3 keV. This component 
is best described by a very flat, T ~ power law, 
accounting for 20% of the total 0.5-10 keV flux. The 
hard component is also independently confirmed by 
the Chandra exposure. 

— The physical origin of the hard X-ray tail in the halo is 
not clear. We argue that contribution from a popula- 
tion of faint Low Mass X-ray Binaries in unlikely. If due 
to IC/CMB, the weak diffuse radio emission and the X- 
ray flux imply a very low magnetic field, << 1/iGauss, 
or a very narrow energy distribution of the high-energy 
electrons. Non-thermal bremsstrahlung from subrela- 
tivistic particles implies the presence of a large number 
of energetic electrons. 

— The jet X-ray spectrum is described by a power law 
with photon index T ~ 1.2. The thermal model advo- 
cated on the basis of previous ROSAT data (Mack et 
al. 1997a) is formally ruled out by the EPIC and ACIS 
data. However, based on the limited signal-to-noise ra- 
tio of both datasets, the possibility that a fraction of 
the jet X-ray emission is thermal remains open. 

— An IC/CMB interpretation of the jet X-ray emission 
implies strong beaming (5 ~ 10) at equipartition. 
Alternatively, weak beaming is possible if the magnetic 
field in the jet is more than one order of magnitude be- 
low equipartition. A similar conclusion holds from the 
Dermer & Atoyan (2002) model. 



A glance at the Chandra image reveals the presence 
of several point sources in the field of the target. We used 
the CIA0 tool wavdetect to search for serendipitous X-ray 
sources in the field of view. In the algorithm, the param- 
eter scale (a list of radii, in image pixels, of Mexican 
Hat wavelet functions) was left free to range between 1 
and 16, and the parameter threshold (the number of de- 
tected spurious sources in a pixel map) held fixed at 10~ 9 . 
The algorithm returns a list of elliptical regions that define 
the positions and the shapes of the detected sources. We 
next used the coordinates from wavdetect and its associ- 
ated error regions to search for their optical counterparts 
on ESO archival plates. We used circular search regions 
with radii of 5" because it is known that most reprocessed 
ACIS-I observations have an offset of up to 1.5". 3 The 
properties of the point sources are listed in Table 1. 

Table 1. Field Sources 



R.A. 


DEC. 


X 


Y 


Rate (c/s) 


Notes 


248.09961 


82.47877 


4177 


3596 ± 21.02 


4.69 




247.97225 


82.48873 


4299 


3669 ± 13.75 


3.74 


(a) 


247.80649 


82.54228 


4455 


4062 ± 5.83 


2.45 




247.67471 


82.53810 


4580 


4032 ± 34.00 


5.92 




248.36962 


82.53544 


3920 


4011 ± 4.80 


2.24 




248.33852 


82.53210 


3949 


3986 ± 16.65 


4.12 




248.25310 


82.52024 


4030 


3899 ± 5.80 


2.45 




248.23887 


82.52276 


4044 


3918 ± 4.81 


2.24 




248.51997 


82.51474 


3776 


3860 ± 4.87 


2.24 




248.16867 


82.57664 


4111 


4312 ± 7.68 


2.83 





Columns: l=Right Ascension at J2000; 2=Declination at 
J2000; 3, 4=Physical Coordinates on the detector; 5=Count 
rate and error in the energy range 0.3-9 keV. 

(a) optically identified, U1650-01867603, RA=247.9735, 
DEC=+82.48894, i?=17.5 mag, B=19.0 mag, offset from 
X-ray source =0.016'. 

(b) optically identified, U1725-00533708 (NGC 6252), 
RA=248.16904, DEC=+82. 57678, R=9.2 mag, B=11.4 mag, 
offset from X-ray source =0.008'. 

Outside the extraction radius of 4' there are two more optically 
identified sources: 

(cl) U1650-01870981 , RA=248.63446, DEC=+82.48589, 
i?=12.1 mag, B=12.6 mag; 

(c2) U1725-00532074 (2MASS J16290730+8231573), 
RA=247.27754, DEC=+82.53217, i?=15.5 mag, B=16.6 
mag. 
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Fig. 3. Top: EPIC pn+MOS spectra of the halo of 
NGC 6251, extracted in an annulus of radii 35-215". In 
the top panel, the top spectrum is the pn, the bottom two 
datasets are the MOS1 and MOS2 data. The data were fit- 
ted with a thermal model; excess flux is apparent above 3 
keV, suggesting a harder component. Bottom: Confidence 
contours at 68%, 90%, and 99% confidence for the tem- 
perature and abundance of the thermal component, from 
the fits to the EPIC pn+MOS data. 



Radius [arcsec] 

Fig. 2. Radial profiles of NGC 6251. (a, Top): EPIC 
MOS2 profile in total energy band 0.3-10 keV compared 
to the instrumental PSF (solid line; see text) and back- 
ground (dashed line), (b, Middle): ACIS-I radial profile, 
compared to the instrumental PSF (solid line). A clear ex- 
cess over the PSF is present between 20-100". ( c, Bottom): 
Comparison of the EPIC MOS2 profile at soft (0.3-2 keV, 
solid line plus errors) and hard (2-10 keV, open diamonds) 
X-rays. The soft and hard X-ray profiles are similar. 
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Fig. 4. Plot of the predicted synchrotron flux at 326 MHz 
versus magnetic field for the halo, assuming two different 
slopes for the radio spectrum. The horizontal dotted and 
dashed lines are the upper limits to the halo and total 
fluxes at 326 MHz from the WSRT map of Mack et al. 
(1997b). 




15 

Log u(Hz) 

Fig. 5. Spectral energy distribution of the outer jet. The 
dotted line is the extrapolation of the radio spectrum, 
which overestimates the measured X-ray flux by more than 
two orders of magnitude. 



Fig. 6. Plot of the jet Doppler factor, 8, versus the mag- 
netic field, B, as expected from the IC/CMB process 
(dashed line), SSC process (continuous line), and from the 
equipartition condition (dotted line). Assuming equipar- 
tition, an X-ray origin via IC/CMB implies beaming, 
8 ~ 10. Weak beaming is possible provided the magnetic 
field is very small, ~ 0.1/xGauss, and out of equipartition. 
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Fig. 7. Fit to the radio-to-X-ray spectral energy distribu- 
tion of the jet with the Dermer & Atoyan (2002) model. 
We assumed that the Lorentz factor of cooled parti- 
cles (marking the spectral break at 10 11 Hz) is 7 coo i = 
1.5 x 10 5 . Particles are injected at a rate Q = 7 x 10 49 



sec , with slope p — 2.28. 



